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Abstract 

Previous generations of X-ray observatories revealed a group of massive binaries that were relatively bright X-ray emitters. 
This was attributed to emission of shock-heated plasma in the wind-wind interaction zone located between the stars. 
With the advent of the current generation of X-ray observatories, the phenomenon could be studied in much more detail. 
In this review, we highlight the progress that has been achieved in our understanding of the phenomenon over the last 
15 years, both on theoretical and observational grounds. All these studies have paved the way for future investigations 
using the next generation of X-ray satellites that will provide crucial information on the X-ray emission formed in the 
innermost part of the wind-wind interaction. 
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1. Introduction 

Early-type stars of spectral type OB or Wolf-Rayet 
(WR) feature powerful, highly supersonic stellar winds 
that are driven by the radiation pressure in numerous spec¬ 
tral lines. These winds combine large mass-loss rates and 
wind velocities and thus they vehicle important quantities 
of kinetic power. When two such stars are bound together 
by gravity in a binary system, their winds interact and part 
of the kinetic energy is converted into heat. This interac¬ 
tion gives rise to a number of observational signatures that 
span a wide range of the electromagnetic spectrum, from 
radio waves to the 7-ray domain. In this review, we focus 
on one of the most spectacular observational consequences 
of colliding winds: the X-ray emission that is produced by 
the shock-heated plasma in the wind interaction zone. 

Two years before the first detection of X-ray emis¬ 
sion from massive stars, Prilutskii and Usov (1976) and 
Cherepashchuk (1976) proposed that collisions between 
the wind of one component of a massive binary system 
and either of the wind, photosphere or magnetosphere of 
the companion star should produce strong X-ray emission. 
Observations with EINSTEIN, and lateron with ROSAT, 
revealed that single and binary early-type stars are mod¬ 
erate X-ray emitters. However, for the binary systems, the 
observed flux level turned out to be significantly lower than 
expected from early theoretical considerations. Still, it was 
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found that many of the brighter X-ray sources among mas¬ 
sive stars were binary systems, and that, on average, mas¬ 
sive binary systems had a larger Lx/Uboi ratio than single 
massive stars. This conclusion held for binaries hosting 
Wolf-Rayet stars (Pollock, 1987) as well as for systems of 
spectral-type 0-1-0 (Chlebowski and Garmany, 1991). 
Moreover, the level of X-ray emission from some of these 
systems was found to be variable (Corcoran, 1996). These 
observational results were interpreted in terms of the col¬ 
lision of the stellar winds and led to the general paradigm 
that colliding wind binaries should be X-ray brighter than 
single massive stars. 

This rather simple picture has been refined quite sub¬ 
stantially over the last two decades as more detailed obser¬ 
vations and more sophisticated models became available. 
In the following we will summarize the progress that has 
been achieved in understanding the phenomenon both the¬ 
oretically (Sect. 2) and from an observational point of view 
(Sect. 3). We deliberately leave 77 Car aside as this system 
is subject of a dedicated review in this volume (Corcoran 
et ah). In Sect. 4, we highlight our conclusion and discuss 
some future perspectives. 

2. Theoretical models of wind-wind interactions 

Following the observational detection of X-ray emis¬ 
sion from colliding winds at a flux level significantly be¬ 
low the anticipated one, major efforts have been devoted 
to improve theoretical models of wind-wind interactions. 
Whilst some analytical or semi-analytical formalisms have 
been proposed, most of the work has focused on the devel¬ 
opment of hydrodynamic simulations. Over recent years. 
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the progresses in numerical techniques and computational 
resources allowed for numerical simulations of growing so¬ 
phistication that account for an increasing number of phys¬ 
ical processes. In this section, we first present some general 
considerations, before we discuss the more recent develop¬ 
ments that resulted in the current theoretical picture of 
wind-wind interactions. 


2.1. The basic picture 

The wind interaction zone is contained between two 
oppositely faced shocks separated by a contact discontinu¬ 
ity. The location and shape of this contact discontinuity 
are mainly determined by the value of the wind momen¬ 
tum ratio (e.g. Shore and Brown, 1988; Luo et ah, 1990; 
Stevens et ah, 1992; Usov, 1992; Canto et ah, 1996): 


V = 


Ml vi 
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where Vj and Mj stand for the pre-shock velocity and 
mass-loss rate of star j. To first approximation, the in¬ 
ner part of the wind interaction zone can be described as 
a shock cone wrapped around the star with the weaker 
wind. A popular method to compute the shape of the 
contact discontinuity is the analytical formalism proposed 
by Canto et al. (1996). These authors study the prop¬ 
erties of the wind interaction zone for winds colliding at 
their terminal velocities assuming that the wind interac¬ 
tion zone forms a thin shell around the contact disconti¬ 
nuity. Using first principles. Canto et al. (1996) provide 
non-linear equations for the shape of the contact disconti¬ 
nuity, as well as expressions for the mass surface density 
of the wind interaction zone and the tangential velocity 
along the contact discontinuity. 


At the shock front, the kinetic energy of the inflow¬ 
ing wind normal to the shock is transformed into thermal 
energy. Given the typical velocities of the stellar winds of 
early-type stars, the shocked material is heated to temper¬ 
atures of several 10^ K. What happens to the shock-heated 
plasma depends strongly on the efficiency of radiative cool¬ 
ing. Stevens et al. (1992) propose to quantify the efficiency 
of this process via the parameter 
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where tcooi and tesc stand for the cooling and escape time, 
respectively^. Here v, d and M indicate the pre-shock 
velocity (in 1000kms“^), the separation between the star 
and the shock (in 10^ km) and the wind mass-loss rate (in 
10“^ Mq yr“^), respectively. Each of the two winds has its 
own value of x (Stevens et al., 1992). For situations where 
X << 1, the cooling is very fast and the shock-heated 


^Equation 2 assumes that the cooling function is independent of 
the temperature which is only the case over a limited temperature 
range near 10^ K, see e.g. Fig. 10 of Stevens et al. (1992). 


material instantaneously radiates its energy away so that 
the wind is essentially isothermal. Conversely, if x > 1; 
the shocked plasma is adiabatic. In this case, radiative 
cooling is absent and the shock-heated plasma in the wind 
interaction zone only cools through adiabatic expansion. 

Luo et al. (1990) presented the first hydrodynamic sim¬ 
ulation of colliding winds. These authors focused on two 
extreme situations: a purely radiative and a purely adia¬ 
batic system. The hydrodynamic simulations of Stevens 
et al. (1992) provided a major step forward as they were 
the first to explicitly include radiative cooling in a self- 
consistent way, therefore allowing to consistently treat sys¬ 
tems that are partially radiative. 

An important limitation of these early simulations con¬ 
cerns their spatial resolution in the treatment of radiative 
wind interaction zones. Indeed, whilst adiabatic shocks are 
self-similar and simply scale with the separation between 
the stars, the situation is more complex when radiative 
cooling is important. In the latter case, the cooling length 
introduces another scale length in the problem (Stevens 
et al., 1992). Dealing with a spatial grid that fully resolves 
the structures of radiative shocks whilst simultaneously en¬ 
compassing the whole wind interaction zone is very chal¬ 
lenging in terms of numerical resources. Early hydrody¬ 
namic simulations thus failed to resolve the cooling layers. 
Nowadays the situation has improved thanks to the use of 
adaptive mesh refinement algorithms (e.g. Parkin et al., 
2014, and references therein), although highly radiative 
shocks remain a major challenge for numerical simulations. 
These difficulties of the hydrodynamic simulations to cope 
with radiative shocks motivated the development of a semi- 
analytic formalism for computing the X-ray emission from 
such shocks (Antokhin et al., 2004). In this approach, 
the wind acceleration is treated via the wind acceleration 
laws of single stars. The wind interaction is idealized as 
an axisymmetric, laminar, steady radiative shock and the 
cooling layer behind the shock is assumed to be isobaric. 
Since mixing of hot and cool material in the post-shock 
region is ignored in this model, it provides an upper limit 
on the actual level of X-ray emission. 

An important question is indeed how much X-ray emis¬ 
sion is produced by the shock-heated plasma. Usov (1992) 
presented a number of analytical recipes to predict the X- 
ray luminosities of winds colliding at Voo ■ When compared 
to actual observations, this formalism was frequently found 
to overpredict the level of the X-ray flux. Stevens et al. 
(1992) showed that in the adiabatic case, one expects the 
X-ray luminosity to approximately scale as 
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where the symbols have the same meaning as before. Con¬ 
versely, for the purely radiative case, Stevens et al. (1992) 
found that 
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where / is the fraction of the wind involved in the collision. 

The question of which shocked wind is responsible for 
most of the X-ray emission in colliding wind interactions 
was re-addressed by Pittard and Stevens (2002). They 
showed that in situations where both winds cool rapidly, 
the faster wind is generally the one that dominates the X- 
ray emission. Conversely, in adiabatic systems, the stronger 
wind (i.e. the one with the larger wind momentum, M v) is 
usually the one that emits most of the X-rays. The conclu¬ 
sions of Pittard and Stevens (2002) largely disagree with 
the recipes of Usov (1992). In general, both the analytic 
and hydrodynamic calculations from these early models 
predict X-ray luminosities that are an order of magnitude 
too large compared to the observed values. Stevens et al. 
(1992) also pointed out that their models for V 444 Cyg 
(WN5 -I- 06, e = 0, Porb = 4.2 d) predict a plasma tem¬ 
perature well above that measured on EINSTEIN spectra. 

2.2. The physics of the shoeks 

Our current picture of hydrodynamic shocks in collid¬ 
ing wind binaries is significantly more complex than a sim¬ 
ple application of the Rankinc-Hugoniot conditions in the 
limit of strong shocks. Additional effects need to be taken 
into account. These concern modifications of the proper¬ 
ties of the pre-shock flow as well as microscopic effects in 
the post-shock plasma. 

In a massive binary system, the stellar wind material 
interacts with the radiation fields and the gravity of both 
stars. As a result, applying the stellar wind formalism of 
single massive stars to the winds in binary systems only 
provides a first order approximation. Stevens and Pollock 
(1994) were the first to investigate this problem along the 
axis between the stars in close massive binary systems. 
They found that the radiation field of the companion in¬ 
hibits the acceleration of one star’s wind preventing it from 
reaching Voo- Whilst the gravitational field of the com¬ 
panion compensates part of the effect, it cannot prevent 
a lowering of the pre-shock velocity. Radiative inhibition 
therefore lowers the post-shock temperature and Stevens 
and Pollock (1994) suggest that this effect could explain 
why the X-ray emission from close massive binaries is sig¬ 
nificantly softer than expected from the terminal velocities 
of the stellar winds. 

Gayley et al. (1997) considered another consequence 
of the radiation pressure of the companion. This effect, 
called radiative braking, is most relevant for systems that 
exhibit a large imbalance between the wind momenta, such 
as WR -I- O binaries. In such systems, the dominant wind 
has reached an asymptotic velocity, but is suddenly de¬ 
celerated by the radiation it encounters upon approach¬ 
ing the companion. Radiative braking could lead to a 
wider opening angle of the shock cone (e.g. in V 444 Cyg, 
Sect. 3.2.1). The same effect could prevent a collision of 
the dominant wind onto the companion’s surface in situa¬ 
tions where the companion’s wind is unable to sustain the 


momentum of the incoming wind (e.g. at periastron pas¬ 
sage in 72 Vel, WC8 -k 07.5III/9I, e = 0.33, Porb = 79d). 
Gayley et al. (1997) provide a semi-analytical recipe to es¬ 
timate whether or not radiative braking is important for 
a given system. 

Another way the radiative acceleration of the pre-shock 
wind could be altered was studied by Parkin and Sim 
(2013). These authors introduce the concept of self-regulated 
shocks: the X-ray radiation from the shocked material ion¬ 
izes the wind material ahead of the shock, thereby sup¬ 
pressing the line force and inhibiting the wind acceler¬ 
ation. Parkin and Sim (2013) find that this effect can 
lead to a lower pre-shock velocity, a correspondingly lower 
temperature of the post-shock plasma, and enhanced cool¬ 
ing, possibly also rendering the shock front unstable. Yet 
the intrinsic X-ray luminosity of the shocks was found to 
remain largely unaltered. For systems with very differ¬ 
ent wind momenta, Parkin and Sim (2013) find that this 
effect renders radiative braking less effective and signifi¬ 
cantly extends the range of binary separations over which 
a wind-photosphere collision is likely to occur. 

Myasnikov and Zhekov (1998) considered the effect of 
thermal conduction due to electrons on the properties of 
the shocks in an adiabatic wind interaction zone. Accord¬ 
ing to their results, thermal conduction leads to preheat¬ 
ing zones in front of the shocks where the temperature is 
higher than expected. This in turn leads to lower Mach 
numbers of the inflowing gas. Myasnikov and Zhekov (1998) 
concluded that the higher the efficiency of thermal con¬ 
duction, the lower the temperature of the shocked gas be¬ 
cause the energy is distributed over a larger quantity of 
material. They predict that the resulting X-ray spectrum 
should be softer than in the absence of conduction. This 
thermal conduction is not to be confused with numerical 
heat conduction in hydrodynamical simulations. Parkin 
and Pittard (2010) showed that the dynamics of the hot 
gas and the corresponding X-ray emission are affected by 
this purely numerical problem to an extent that depends 
on the density and temperature contrast across the con¬ 
tact discontinuity. The strongest effects were found when 
the shocked wind of one star behaves quasi-adiabatically 
whilst the shocked wind of the other star is strongly ra¬ 
diative. Higher resolution simulations reduce the artificial 
heat transfer, but increasing the resolution is not always 
practicable. 

Energy losses via Gomptonization, i.e. the redistribu¬ 
tion of energy from the thermal electrons in the shocked 
plasma to the photons of the stellar radiation field follow¬ 
ing inverse Gompton scattering, also affect the properties 
of the wind interaction region. Myasnikov and Zhekov 
(1991, 1993) showed that energy losses through Gomp¬ 
tonization lead to an increase of the density and a de¬ 
crease of the temperature of the gas in the interaction zone. 
Therefore, whilst the emission measure of the shocked gas 
increases, the peak energy of its emission shifts towards 
lower energies. As a result. X-ray emission gets softer and 
thus more sensitive to absorption by the cool wind mate- 
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rial. 

Finally, in colliding wind binary systems where rela¬ 
tivistic particles are accelerated through diffusive shock 
acceleration (see Sect. 2.6), the pre-shock velocity can also 
be modified through so-called shock modification (Pittard 
and Dougherty, 2006). This effect results from the diffu¬ 
sion of non-thermal ions upstream of the shock. These rel¬ 
ativistic ions then exert a back-pressure on the pre-shock 
flow leading to a decrease of the pre-shock velocity and 
the formation of a so-called shock precursor (Pittard and 
Dougherty, 2006). The reduced velocity jump leads to a 
softer X-ray spectrum. In eccentric binary systems, the 
degree of shock modification could be dependent on the 
orbital phase. 

In the strong shock limit, the kinetic energy normal to 
the shock is converted into heat according to (e.g. Stevens 
et ah, 1992) 

kTs = (5) 

where Tg, nip and v±^ are the post-shock temperature, the 
particle mass and the particle’s pre-shock velocity perpen¬ 
dicular to the shock, respectively. Because of the very 
different particle masses, electrons and ions are thus ex¬ 
pected to display quite different post-shock temperatures. 
As pointed out by several authors (Usov, 1992; Zhekov and 
Skinner, 2000; Zhekov, 2007; Pollock et ah, 2005, amongst 
others), in wide binary systems, the postshock densities 
are such that the timescale for ion and electron tempera¬ 
ture equalization through Coulomb interactions are com¬ 
parable to the escape times. Therefore, one would expect 
an electron temperature much lower than the ion tempera¬ 
ture. Yet, the observed values of the electron temperature 
are not that low, but rather fall in the range between 0.1 
and 1.0 times the ion temperature^. This led to the con¬ 
clusion that the shocks must be collisionless (Zhekov and 
Skinner, 2000; Zhekov, 2007; Pollock et ah, 2005). How¬ 
ever, the understanding of this phenomenon in the context 
of colliding wind binaries is still fragmentary, and no clear 
rule how to deal with this effect has emerged so far. For in¬ 
stance, Pittard and Dougherty (2006) suggest that shocks 
that are efficient in accelerating particles to relativistic en¬ 
ergies might equilibrate the electron and ion temperatures 
much faster than shocks where no such acceleration takes 
place. 

In a similar context, many observational studies and 
early models assumed that the shock-heated plasma is 
in ionization equilibrium. However, in a some cases the 
timescales needed to reach ionization equilibrium through 
collisions are longer than the dynamical timescales and the 
plasma is thus likely not in equilibrium. Zhekov (2007) 


^As an example, Pittard (2009) assumed in his models that the 
electron temperature immediately behind the shock amounts to 0.2 
times the mean plasma temperature. 


proposed to quantify the importance of non-equilibrium 
ionization effects via the criterion 

rNEi = i.2i|^ (6) 

Here, Xe is the ratio between the electron and ion number 
densities, JI is the mean atomic weight for ions, d is the 
distance from the star to the shock in 10^^ km, M is the 
mass-loss rate in 10 “® Mq yr“^ and v the wind velocity in 
units 1000 km s“^. For Fnei ^ 1, non-equilibrium effects 
should play a minor role, whereas they are expected to be 
important if Fnei < 1 (Zhekov, 2007). 

2.3. The impact of orbital motion 

Since the first 2-D models of colliding winds which ne¬ 
glected the orbital motion and assumed an axisymmetric 
wind interaction region, there has been enormous progress 
in the treatment of the impact of orbital motion on the 
wind-wind interaction. Currently, state-of-the art models 
use 3-D adaptive mesh refinement hydrodynamical codes 
accounting for radiative driving, gravity, radiative cooling 
and orbital motion (e.g. Parkin et ah, 2011b, 2014). 

Pioneering work on this subject was performed by the 
Zurich group (e.g. Walder, 1995; Walder et ah, 1999; Folini 
and Walder, 2000) and by Pittard and Stevens (1999). For 
instance, Walder et al. (1999) presented a 3-D hydrody¬ 
namic simulation of 7 ^ Vel performed with the goal to 
reproduce the observed AOS’AT lightcurve of the system 
(Willis et ah, 1995). Although the models of Walder et al. 
(1999) assumed a purely adiabatic interaction, whereas the 
real one is at least partially radiative, the agreement be¬ 
tween the simulated and observed lightcurves was rather 
satisfactory. In these models the unperturbed winds were 
assumed to have a constant velocity equal to Voc ■ This as¬ 
sumption was still used in more recent work (e.g. Lemaster 
et ah, 2007). 

In the same context, Okazaki et al. (2008) performed 3- 
D smoothed particle hydrodynamics simulations to model 
the RXTE X-ray lightcurve of 7 Car (Porb = 5.54yr, e ^ 
0.9). Okazaki et al. (2008) assumed that the bulk of the 
X-ray emission arises at the apex of the shock cone. Their 
models then provided the density structure of the wind in¬ 
teraction zone, hence allowing to compute the associated 
variations of the column density along the line-of-sight to¬ 
wards the apex. 

At the same time, Parkin and Pittard (2008) presented 
a simple 3-D dynamical model to evaluate the orbital- 
induced curvature of the colliding wind region. In this 
work, the winds were again assumed to collide at terminal 
velocity. The apex of the colliding wind region was com¬ 
puted adopting the formalism of Canto et al. (1996), whilst 
the downstream contact discontinuity was obtained assum¬ 
ing a ballistic flow. The aberration of the apex due to the 
orbital motion was approximated by tan p, = (with Vi^ 
the slower of the two wind velocities). Parkin and Pittard 
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(2008) used this model mainly to estimate the X-ray at¬ 
tenuation by the circumstellar material. The same model 
was subsequently applied by Parkin et al. (2009) to model 
the RXTE lightcurve of t] Car. In contrast with the work 
of Okazaki et al. (2008), Parkin et al. (2009) accounted for 
the spatial extension of the X-ray emission zone. Still, the 
model failed to reproduce the broad minimum observed 
in the lightcurve. Parkin et al. (2009) accordingly sug¬ 
gested that the secondary wind suffers radiative inhibition 
at periastron, implying a lower wind velocity and hence 
stronger post-shock cooling. As a result, the wind colli¬ 
sion zone could collapse onto the secondary (see also the 
review by Corcoran et al., this volume). 

Pittard (2009) presented 3-D hydrodynamical models, 
incorporating gravity, the driving of the winds (including 
radiative inhibition), the orbital motion of the stars and 
radiative cooling of the shocked plasma. He applied these 
models to a series of O -I- O case studies. These sim¬ 
ulations revealed some interesting features. For instance, 
Pittard (2009) showed that, in systems with equally strong 
winds {r] = 1), the gas near the leading shocks is less dense 
and cooler than the gas near the trailing shocks. This is 
because, at a given downstream distance from the binary 
axis, the leading shock is more oblique than the trailing 
shock. 

Finally, Lamberts et al. (2012) considered the structure 
of an adiabatic colliding wind region on scales much larger 
than the orbital separation. From first principles, the or¬ 
bital motion is expected to turn the shock structure into 
a spiral. Lamberts et al. (2012) used an adaptive mesh 
refinement technique to simulate the large scales whilst si¬ 
multaneously keeping a high resolution close to the binary. 
They found that Kelvin-Helmholtz instabilities are trig¬ 
gered by the effect of the orbital motion on the adiabatic 
wind interaction zone even for winds of equal strength. In 
some cases, i.e. when there exists a steep velocity gradient 
between the winds, the growth of this instabilities prevents 
the formation of a stable spiral structure (Lamberts et al., 
2012 ). 

2.4- The ejfects of instabilities and clumps 

Stevens et al. (1992) drew attention to the importance 
of dynamical instabilities in their hydrodynamic simula¬ 
tions. They reported the presence of Kelvin-Helmholtz in¬ 
stabilities which arise from the shear when the velocities of 
the colliding winds are different (see also Lamberts et al., 
2011). Stevens et al. (1992) also reported on thin-shell 
instabilities which form when radiative cooling produces 
shells of cold dense gas behind the shock. 

An important question is of course whether or not these 
instabilities and their growth reflect a physical phenomenon 
or stem from limitations of the numerical scheme or spe¬ 
cific model assumptions. This aspect was addressed by 
Myasnikov et al. (1998) for systems with highly radia¬ 
tive wind interaction zones. In simulations of such sys¬ 
tems, its is frequently assumed that the energy provided 
by stellar radiation prevents the shocked material from 


having its temperature drop below some threshold value. 
However, this threshold value is usually not derived self- 
consistently. Myasnikov et al. (1998) therefore performed 
numerical simulations with three different values of this 
threshold temperature: 10®, 10® and 10^ K. These authors 
found that the choice of the adopted value had a strong 
impact on the development of the instabilities in the wind- 
wind interaction zone. 

In radiative wind-wind interactions, the dynamical in¬ 
stabilities determine the shape and size of the wind interac¬ 
tion zone as well as its small-scale structure (Walder and 
Folini, 2000). Stevens et al. (1992) predicted that these 
turbulent motions could lead to variability of the X-ray 
flux of order 10%. But these instabilities could also have 
other far-reaching consequences. Indeed, Kee et al. (2014) 
noted that the thin-shell instabilities transform the com¬ 
pressive shock into highly oblique shocks, thus reducing 
the heating of the material and therefore the associated X- 
ray emission. These authors presented a parameter study 
of such instabilities showing that in the radiative regime, 
the X-ray emission follows indeed a linear scaling with 
mass-loss, as expected from Equation (4), but reduced by 
about a factor 1/50 compared to analytic radiative-shock 
models without instabilities (similar to the formalism of 
Antokhin et al., 2004). Though the situations studied by 
Kee et al. (2014) consistently yield this factor 1/50 reduc¬ 
tion, the universality of this reduction factor and/or its 
dependence on key parameters such as the Mach number 
of the pre-shock flow still need to be established. 

There is nowadays strong observational evidence that 
the stellar winds of massive stars are highly clumped (e.g. 
Eversberg et al., 1998). Cherepashchuk (1990) tried to ex¬ 
plain the low luminosity of colliding wind systems com¬ 
pared to theoretical models, by the fact that only the 
smooth inter-clump wind is affected by the collision zone 
between the stars, the clumps passing unaffected through 
it (and later colliding onto the companion, without X- 
ray production). Aleksandrova and Bychkov (2000) ar¬ 
gued that the shocked clumps should be highly radiative, 
thereby counterbalancing the 1/d trend expected for the 
adiabatic behaviour of the lower density smooth wind com¬ 
ponent in eccentric systems such as WR 140 (WC7 -I- 05.5, 
e = 0.89, Porb = 7.93yr). However, these conclusions 
are not supported by the hydrodynamic simulations of the 
collision between clumpy winds in WR 140 presented by 
Pittard (2007). In the latter work, it was shown that, in 
situations where the wind interaction zone is adiabatic^, 
the clumps are rapidly destroyed in the post-shock region. 
This clump destruction occurs mainly through dynamical 
instabilities such as the Kelvin-Helmholtz instability. Pit¬ 
tard (2007) accordingly concluded that the X-ray emission 
of the wind interaction zone is insensitive to clumping. 
Still, the clumps penetrating the wind interaction zone 


^So far, no study of the influence of clumps on radiative wind 
interaction zones has been presented. 
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Figure 1: 3-D hydrodynamic simulation of the wind-wind interaction in WR 22 including the effect of radiative acceleration (from Parkin 
and Cosset, 2011). Left: gas density in the orbital plane at three different orbital phases. Right: temperature of the gas in the orbital plane. 
The WR star is the star with the stronger wind. Each plot covers a region of 1.2 x 10^“^ cm on one side. Credit: Parkin, E.R. , Sz Cosset, E., 
ASzA, 530, A119, 2011, reproduced with permission ©ESO. 
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trigger turbulent motion that slightly puffs up the wind 
interaction zone. 

2.5. Orbital variability 

Using the same code as Stevens et al. (1992) (i.e. as¬ 
suming axisymmetric winds colliding at their terminal ve¬ 
locity), Pittard and Stevens (1997) computed synthetic 
X-ray lightcurves and spectra for three different values of 
the wind momentum ratio. They considered two mostly 
adiabatic systems and one radiative case. By construction, 
these lightcurves were symmetric about phase 0.5. Most 
of the predicted variations were hence due to the chang¬ 
ing optical depth along the line of sight (i.e. larger optical 
depth when the star with the denser wind is in front, de¬ 
creasing the observed X-ray flux) though occultation by 
the stellar bodies also played a role. 

These results were revised by Pittard and Parkin (2010) 
who used the 3-D hydrodynamical models of Pittard (2009) 
to predict the variations of the X-ray emission from the 
wind interaction regions for different cases of O -|- O bina¬ 
ries. In accordance with Pittard and Stevens (1997), sys¬ 
tems with circular orbits were predicted to exhibit phase- 
dependent X-ray variability only due to changes of the 
line-of-sight column densities. Since regions emitting X- 
rays of different energies are predicted to have different 
spatial extents, the column densities depend on the en¬ 
ergy. Overall, in such systems, the maximum X-ray flux is 
expected at quadrature phases. However, in contrast with 
the models of Pittard and Stevens (1997), the lightcurves 
of short-period circular systems as computed by Pittard 
and Parkin (2010) are not fully symmetric due to orbital 
aberration. Moreover, in these short-period systems ra¬ 
diative cooling is quite efficient and the models hence pre¬ 
dict sharp dips in the soft-band lightcurve when the dense 
thin layer of cooled post-shock gas crosses the line-of-sight. 
The latter feature disappears for longer-period systems 
with adiabatic wind interactions, for which the synthetic 
lightcurves become more symmetric. 

For eccentric systems, Pittard and Parkin (2010) pre¬ 
dict also changes in the intrinsic emission. However, in 
the systems considered by Pittard and Parkin (2010), the 
behaviour is far more complex than the 1/d relation ex¬ 
pected from Equation (3). This is because the wind colli¬ 
sion region changes from adiabatic at apastron, to highly 
radiative at periastron. Moreover, the pre-shock veloci¬ 
ties also change. As a result, Pittard and Parkin (2010) 
predict that the soft emission reaches a maximum at pe¬ 
riastron (</) = 0.0), whilst the hard emission peaks before, 
and then drops to almost its minimum value at periastron. 
Moreover, in the 3-D models of Pittard (2009) some dense 
clumps of cool material formed during periastron passage 
were found to survive inside the wind interaction region. 
Because of their high density and inertia, it takes some 
fraction of the orbital cycle (about half the orbital period 
in the simulations of Pittard, 2009) to push them out. As 
a result of the history of the wind interaction zone, the 
results of Pittard and Parkin (2010) thus reveal strong 


hysteresis effects with the emission being softer after pe¬ 
riastron. These effects are reminiscent of the hysteresis 
found in the synthetic thermal radio lightcurves computed 
by Pittard (2010). 

Even more complex phenomena have been predicted for 
eccentric systems as a result of instabilities of highly radia¬ 
tive shocks. In their effort to explain the X-ray lightcurve 
of WR 22 (WN7 -f O, e ~ 0.56, Porb = 80 d). Parkin and 
Cosset (2011) presented hydrodynamic simulations indi¬ 
cating that the collision region is very unstable if the ac¬ 
celeration of the winds is taken into account. Parkin and 
Cosset (2011) show that, as the stars approach perias¬ 
tron, the ram pressure of the WR wind overwhelms that 
of the 0-star and following a disruption of the shocks by 
non-linear thin-shell instabilities, the wind collision region 
collapses onto the O star. 

Finally, synthetic X-ray line profiles formed in the hot 
post-shock gas and their orbital variations were computed 
by Henley et al. (2003). These computations were based 
on the density and velocity fields from snapshots of 2-D 
axisymmetric hydrodynamic simulations of adiabatic wind 
collisions assuming that the winds collide at their terminal 
velocity. Henley et al. (2003) focused on profiles of the Lya 
transitions of O, Ne, Mg, Si and S which fall into the sen¬ 
sitivity range of the current generation of high-resolution 
X-ray spectrographs. However these lines are also emitted 
by the hot plasma intrinsic to the winds of the individual 
stars that make-up the binary. Therefore, their profiles 
are likely contaminated by the intrinsic emission of these 
stars. This is why Rauw et al. (2016) recently computed a 
series of synthetic line profiles for the Fe xxv K line near 
6.7keV which is much less affected by emission intrinsic 
to the stars. In their work, Rauw et al. (2016) use the 
semi-analytical formalism of Canto et al. (1996) to predict 
the complicated morphology of the Fe K line complex for 
different systems and different orbital phases. 

2.6. Magnetic fields and particle acceleration 

Although magnetic fields have become an important 
topic in massive stars research, understanding their im¬ 
pact on wind-wind interactions is a subject that is still 
in its infancy. In the context of colliding winds, mag¬ 
netic fields are notably responsible for the synchrotron 
radio emission observed in some of these systems (e.g. 
Pittard and Dougherty, 2006). However, the requested 
field strengths are actually relatively low. For instance, in 
the case of CygOB2 #9 (05-5.51 -f 03-4 HI, e = 0.71, 
Forb = 2.35yr). Parkin et al. (2014) estimated that the 
surface magnetic field required to account for the observed 
synchrotron radio emission is about 0.3 - 52 G. Such values 
are below the detection capabilities of current spectropo- 
larimetric facilities. The lack of detections of stellar mag¬ 
netic fields in particle-accelerating colliding wind binaries 
(Naze et ah, 2012a; Neiner et ah, 2015) could also indicate 
that the fields responsible for the synchrotron emission 
forms in situ, i.e. inside the wind interaction zone. 
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If we consider that a weak stellar magnetic field exists, 
what consequences are expected for the wind-wind inter¬ 
action? Falceta-Gongalves and Abraham (2012) present 
the hrst full magnetohydrodynamic numerical simulation 
of the wind-wind interaction in 77 Carinae assuming both 
stars have a weak dipolar magnetic field with a polar field 
strength of 1 G. They find that the magnetic field reduces 
the strength of the shock since part of the upstream ki¬ 
netic energy is converted into magnetic rather than ther¬ 
mal energy. Behind the shocks, the density of the magnetic 
field line increases. In an adiabatic case with very large 
Mach numbers, the magnetic held component parallel to 
the shock increases by a factor 4. In a radiative shock, the 
effect is much larger. Yet Falceta-Gongalves and Abraham 
( 2012 ) show that the magnetic held evolution at the shock 
is not linearly correlated to density or thermal pressure. 
Finally, depending on its orientation, a magnetic held can 
reduce quite signihcantly the efficiency of electron conduc¬ 
tion (Orlando et ah, 2008). 

Particle acceleration and the subsequent energy losses 
via synchrotron radiation and inverse Compton scattering 
are also important ingredients in the physics of colliding 
winds. Beside the issue of shock modihcation (Pittard and 
Dougherty, 2006, see Sect. 2.2), the substantial amount of 
energy required to accelerate particles affects the proper¬ 
ties of the shocks and the level of thermal X-ray emission. 
Pittard and Dougherty (2006) present a model of the syn¬ 
chrotron radio emission, the thermal X-ray emission and 
non-thermal (inverse-Compton) X-ray and 7 -ray emission 
of WR 140. They argue that the efficiency of diffusive 
shock acceleration varies with orbital phase and that this 
could explain why the X-ray emission of this system devi¬ 
ates from the l/d relation (see Sect. 3.2.2 below). 

Reitberger et al. (2014) solved the hydrodynamic equa¬ 
tions along with transport equations of energetic parti¬ 
cles accounting for diffusive shock acceleration, radiative 
losses, and adiabatic cooling. In their simulations, none of 
the electron energy spectra reaches energies above 1 TeV 
whilst protons reach energies of several TeV. Orbital mo¬ 
tion, although included at a modest level in this simu¬ 
lation, changes the properties of the leading and trailing 
arms. However, although this is not the subject of the 
present review, we note that, apart from 77 Gar (see Gor- 
coran et al. this volume), to date no colliding wind system 
has been detected at 7 -ray energies. Indeed, Werner et al. 
(2013) investigated data from 24 months of Fermi Large 
Area Telescope (LAT) observations and failed to detect 
emission from any out of seven WR -I- O colliding wind 
binaries ( 7 ^ Vel, WR 70, 125, 137, 140, 146, 147) though 
they were expected to be good candidates. For two of the 
stars (WR 140 and 146), the upper limits are well below 
model predictions of Benaglia and Romero (2003); Pittard 
and Dougherty (2006); Reimer et al. (2006); Reimer and 
Reimer (2007, 2009b, a). 


3. The observational picture 

On the observational side, major progress became pos¬ 
sible thanks to the unprecedented sensitivities and spec¬ 
tral resolution of the Chandra and XMM-Newton missions. 
These observatories allowed monitoring the orbital vari¬ 
ability of the X-ray emission of interacting wind systems. 
Thanks to the longevity of these satellites, this became 
possible also for systems with periods of several years. In 
addition, both spacecraft have been used to perform large 
sensitive surveys of clusters and associations hosting mas¬ 
sive stars, thereby allowing to refine the global picture of 
the X-ray emission of massive stars and massive binaries 
in particular. 

3.1. The X-ray luminosity 

Using data from the ROSAT All-Sky Survey, Berghoe- 
fer et al. (1997) showed that the X-ray luminosity Lx of 
0 -type stars scales roughly as 10 “'^Lboi, thereby confirm¬ 
ing earlier findings with EINSTEIN. However, there was 
considerable scatter around this relation. Moreover, the 
multiplicity of many of the 0 -stars in the sample was not 
known, and for many objects only upper limits on the X- 
ray flux were available. 

Whilst no all-sky survey was performed with neither 
XMM-Newton nor Chandra, both observatories have been 
used to perform surveys of clusters and/or of large mas¬ 
sive stellar populations. These campaigns significantly im¬ 
proved upon the older work. For instance, detections for 
the full 0 -star population in a given cluster were avail¬ 
able, reducing biases. In parallel, optical spectroscopic 
monitoring have tremendously improved the binary cen¬ 
sus, allowing for trends to be better identified. 

It must however be noted that details in processing and 
calibrations as well as assumptions on reddening laws and 
bolometric corrections render direct comparisons of Lx® 
and Lx/Lboi values between different papers sometimes 
difficult (see remarks on a 0.2 dex change in log[Lx/Lboi] 
from such simple effects in Sect. 2.2 of Naze et al., 2011). 

Nevertheless, all homogeneous studies using new, sen¬ 
sitive data led to the same conclusion (Sana et ah, 2006; 
Oskinova, 2005; Antokhin et ah, 2008; Naze, 2009; Naze 


®We note that one must be cautious about the definition of Lx- 
Here, the X-ray luminosity will only be corrected for absorption by 
the neutral interstellar medium. X-ray spectra of massive stars also 
display absorption due to ionized wind material, but its value crit¬ 
ically depends on the signal-to-noise ratio of the data and on the 
assumptions of the models used by the authors, resulting in several 
dex differences on the intrinsic luminosity. Therefore, while they 
could at first be considered as more physical, the meaning of fully 
dereddened luminosities is ambiguous. The well-known degeneracy 
between absorption column and temperature in spectral fits (a hot 
plasma with low absorption providing a fit of equivalent quality to 
a warm plasma with more absorption) further limits the meaning of 
“average” temperatures (for examples, see Naze et al., 2014). Unfor¬ 
tunately, both fully dereddened luminosities and average tempera¬ 
tures are still sometimes the only parameters used (e.g. Gagne et al., 
2012 ). 



et al., 2011; Rauw et al., 2014b): the majority of the known 
O + OB binaries display Lx/L\^oi ratios in the same range 
as supposedly single O stars. This means that not all mas¬ 
sive binaries harbour X-ray bright wind-wind collisions; on 
the contrary, such X-ray bright CWBs appear quite rare, 
whatever their period (Naze, 2009; Naze et ah, 2011). 

Still, on average, binaries have a slightly larger Lx/Tboi 
ratio, by about 0.1-0.2dex, especially in harder energy 
bands (Naze, 2009; Naze et ah, 2011; Rauw et ah, 2014b). 
Although systematic, this difference remains below the 
observed scatter around the average Ax/Aboi relation of 
supposedly single objects. Hence it is not formally signih- 
cant and renders the identification of a binary from its sole 
Ax/Aboi ratio rather difficult. 

This lack of strong overluminosities in binaries, appar¬ 
ently independently of the orbital period or the spectral 
types, remains to be explained. Efficient cooling or radia¬ 
tive braking, as well as reduced mass-loss rates certainly 
play a role. Most probably, X-ray bright colliding wind bi¬ 
naries populate only a small region of the parameter space. 
In addition, in some X-ray bright binary systems, the over¬ 
luminosity is not linked to a wind-wind collision but stems 
from the presence of magnetically confined winds in one of 
the stars (e.g. HD 191612 Naze et al., 2010). 

To make things more complicated, most of the long- 
period massive binaries feature highly eccentric orbits and, 
as a result, the X-ray flux of the colliding winds interaction 
should vary quite significantly between apastron and peri- 
astron. However, as these systems spend most of the time 
far away from periastron, snapshot X-ray observations are 
more likely to catch them in a low emission state. It thus 
appears impossible to draw secure conclusions on binarity 
from a single X-ray snapshot observation. 

The situation is more complex for WR stars for which 
EINSTEIN and ROSAT data did not reveal a well de¬ 
fined Tx/Aboi relation (Wessolowski, 1996). For WR bina¬ 
ries, recent observations have generally confirmed the ear¬ 
lier results and extended them to the Magellanic Clouds 
(MCs) (Guerrero et ah, 2008): binaries appear brighter 
than supposedly single stars. This is particularly the case 
for WC stars, where no single object has been detected in 
the X-ray domain up to now (Oskinova et ah, 2003; Skin¬ 
ner et ah, 2006), with the possible exception of the WC4 
star WR 144 (Rauw et ah, 2014b). For WN binaries, it 
remains generally true that binaries are brighter than sin¬ 
gle stars, although there are exceptions where the emission 
does not seem to come from the wind interaction. For in¬ 
stance, the short-period, eclipsing binary CQ Cep (WN6 + 
09, Porb =1.64 d) displays a constant X-ray emission with 
typical characteristics of single WN objects hence without 
obvious traces of colliding wind emission (Skinner et ah, 
2015). 

A link between X-ray brightness and orbital period was 
proposed by Pollock (1987), with short-period systems be¬ 
ing less bright than those with longer periods, but this has 
not been checked yet using recent data. Zhekov (2012) fur¬ 


ther found that even short-period WR -|- O systems appear 
to be adiabatic, despite expectations of radiative collisions 
based on the high density of the WR wind. This prob¬ 
lem could be solved if the mass-loss is both clumpy and 
reduced compared to older evaluations, a scenario which 
seems very plausible at the present time. Finally, Guerrero 
et al. (2008) suggest that the detected X-ray emission of 
WR stars in the MGs is brighter than that of equivalent 
Galactic objects. If confirmed, this effect could be linked 
to the lower metallicity in the MCs, which leads to a re¬ 
duced absorption. 

Another important point concerns the comparison of 
the X-ray luminosities predicted by theoretical models with 
the observed values. In most studies, the comparison be¬ 
tween theory and observation does not account for the 
X-ray emission intrinsic to the individual stars. In some 
cases, the emission from the wind-wind interaction actu¬ 
ally represents only a minor contribution to the full X-ray 
emission. This is the case of 9 Sgr (03.5V((P)) -I- 05- 
5.5 V((f)), e = 0.70, Porb = 9.1 yr, Rauw et al., 2015). An¬ 
other example is provided by the highly eccentric 09 HI 
-b BlIII binary tOri (e = 0.76, Porb = 29.13 d). Pit- 
tard et al. (2000) analysed two ASCA observations of that 
system taken at periastron and apastron. The X-ray lu¬ 
minosities and spectral distributions were found to be re¬ 
markably similar on both observations, in stark contrast 
with the expected order of magnitude variation for a col¬ 
liding winds interaction. This points towards a colliding- 
wind emission that is unexpectedly weak, and an X-ray 
emission dominated by intrinsic shocks in the winds. 

Even neglecting the contribution of the intrinsic emis¬ 
sion, theoretical models quite frequently overestimate Lx 
by at least an order of magnitude (Zhekov and Skinner, 
2000; Skinner et ah, 2007; Zhekov and Park, 2010b; Zhekov, 
2012; Montes et ah, 2013; Zhekov, 2015). As a result of 
stellar wind clumping, there is nowadays a general concen¬ 
sus that mass-loss rates determined with conventional op¬ 
tical, UV, IR or radio diagnostics might need to be reduced 
by typically a factor ^ 3 (e.g. Puls et ah, 2015). Yet, this 
reduction is generally not sufficient to bring the theoretical 
and observed luminosities of colliding wind systems into 
agreement. In this context, De Becker et al. (2004) anal¬ 
ysed an XMM-Newton observation of HD 159176 (07 V 
-b 07V, e = 0.0, Porb = 3.367d). Hydrodynamic simu¬ 
lations as well as the semi-analytic approach of Antokhin 
et al. (2004) were used to theoretically predict the X-ray 
luminosity due to wind-wind interactions. Adopting wind 
parameters from the literature resulted in a striking dis¬ 
crepancy between the predicted and the observed X-ray 
luminosities. A significant reduction of the mass-loss rates 
of the stars was necessary to alleviate part of the discrep¬ 
ancy. The more recent model of Pittard and Parkin (2010) 
with mass-loss rates of 2 x 10 “^ Mq yr“^ is able to account 
for the emission level of HD 159176, although one should 
stress that this model again does not account for intrinsic 
emission. 
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In the case of CygOB2 ^^^9, Parkin et al. (2014) also 
found that mass-loss rates in their simulations had to be 
reduced by a factor 7.7 to match the observed level of the 
X-ray emission. 

In the case of the trapezium-like system QZ Car which 
hosts two binary systems with 09.71 and 08 III primaries, 
Parkin et al. (2011a) noted from simple analytical consid¬ 
erations that the wind interaction in the binary containing 
the 08 III star should produce an X-ray emission 10 times 
stronger than the one observed with Chandra. Parkin et al. 
( 2011 a) suggest that in this system, the wind interaction 
could be disrupted by mass transfer via Roche lobe over¬ 
flow from the 08 III star towards its companion. 

Another striking example is provided by the study of 
WR 22 by Parkin and Cosset (2011). These authors used 
3-D adaptive-mesh refinement simulations to investigate 
the dependence of the X-ray emission associated with the 
wind-wind collision on the acceleration of the stellar winds, 
radiative cooling, and orbital motion. They accounted for 
radiative driving, gravity, optically-thin radiative cooling, 
and orbital motion. The adopted mass-loss rates were 
1.6 X 10“® and 2.8 x lO“^M 0 yr“^ for the WN7 and the 
09 V star, respectively. Parkin and Cosset (2011) showed 
that when a stable wind collision region exists the mod¬ 
els over-predict the X-ray flux by more than two orders of 
magnitude. They also found that the collapse of the wind 
interaction zone onto the O star substantially reduces the 
discrepancy in the 2-10 keV flux to a factor ^ 6 , though 
the observed spectra were not well matched by the mod¬ 
els. Parkin and Cosset (2011) accordingly suggest that the 
agreement between the models and observations might be 
improved by reducing the mass-loss rates of both stars and 
simultaneously increasing the ratio of the wind momentum 
ratio to a value close to 500 in favour of the WR star. Un¬ 
der these circumstances, no normal wind ram pressure bal¬ 
ance would occur at any orbital phase. In the absence of 
radiative braking this would lead to a permanent collapse 
of the wind collision region onto the O star. It remains to 
be seen whether such drastic changes in the wind parame¬ 
ters are consistent with the observed properties of WR 22 
at other wavelength. 

Beside the overluminosity, another criterion for the pres¬ 
ence of X-ray emission due to colliding winds that is fre¬ 
quently used is the spectral hardness. Indeed, fitting the 
spectra of WR -I- O binaries with X-ray bright wind-wind 
interactions requires usually at least two thermal compo¬ 
nents with kT near 0.6 and 2.-4. keV. The harder com¬ 
ponent is generally subject to higher wind absorption, as 
expected for colliding winds: X-rays of higher energy being 
produced closer to the line-of-centers, hence deeper inside 
the winds. Such a hard component is also sometimes seen 
in presumably single 0 -stars, though with a lower strength 
than in X-ray bright colliding wind systems. More gener¬ 
ally, surveys revealed a harder X-ray emission for massive 
binaries compared to presumably single stars, but the dif¬ 
ference is of limited amplitude (see above). Finally, hard 


X-ray emission is also expected in the case of magnetically 
confined winds (though it is not always observed. Naze 
et ah, 2014, , ud-Doula & Naze, this volume), which may 
blur the picture when searching for colliding wind X-ray 
emission. 

Raassen et al. (2003) pointed out that prominent col¬ 
liding wind systems often display a strong Fe K line near 
6.7 keV in their spectra. This line has thus been used as 
a possible indicator of the presence of wind-wind inter¬ 
actions (e.g. in the case of the massive stars in the core 
of the M17 nebula, Mernier and Rauw, 2013), although 
once more it is also seen in some stars with magnetically 
confined winds. We thus conclude that the presence of a 
hard X-ray emission is not necessarily or uniquely linked 
to binarity. 

3.2. Orbital variations 

One of the best ways to confirm the presence of a col¬ 
liding wind region is through the monitoring of the X-ray 
emission. Indeed, in most cases, wind-wind interactions 
should produce phase-locked variations of the observable 
X-ray flux either because of changing absorption along the 
line-of-sight as the stars revolve around each other and/or 
because of the varying separation in eccentric binaries (see 
Sect. 2.5). The only cases were no variations are expected 
would be for circular orbits seen face-on, which is a statis¬ 
tically rare configuration. 

Variability thus represents a powerful diagnostic of col¬ 
liding winds, but it requires a set of time-constrained ob¬ 
servations, which is demanding in terms of observing time 
and sometimes challenging in terms of scheduling. Never¬ 
theless, the results are worth the effort and several mon¬ 
itoring campaigns have been performed. Before, we re¬ 
view the results that have been obtained, let us stress that 
one can also study variations on other timescales. In par¬ 
ticular, Cherepashchuk (1990) expected “flares” occurring 
when clumps collide but no such short-term events were 
reported. Stochastic variability on rather short timescales 
was also expected as a result of dynamical instabilities in 
radiative wind interaction zones (see Sect. 2.4), but again, 
no significant short-term variability of the X-ray emission 
was found so far (e.g. De Becker et ah, 2004; Sana et al., 
2004; Linder et ah, 2006). 

3.2.1. Orbital variations due to changing absorption and 
eclipses 

Photoelectric absorption of the X-rays emitted by the 
shock-heated plasma occurs primarily in the unshocked 
stellar winds along the line-of-sight towards the observer. 
The associated optical depths depend on the wind prop¬ 
erties as well as the binary characteristics. When the two 
stars feature very different wind properties, such as in WR 
-I- OB systems, absorption effects may lead to particularly 
spectacular variations of the X-ray emission. 

In this context, the prototypical case is 7 ^ Vel. Willis 
et al. (1995) reported a strong increase in the ROSAT 
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count rate around the time when the O star passes in front: 
this was readily explained as an effect of the lower absorp¬ 
tion of the more tenuous 0 -star wind filling the shock cone. 
The width of the peak in the X-ray lightcurve is directly 
related to the cone opening angle, and in this case indi¬ 
cates a quite narrow shock cone (half opening angle of 20 - 
30°). Subsequent observations with ASCA (Stevens et ah, 
1996; Rauw et ah, 2000) and XMM-Newton (Schild et ah, 
2004) revealed that the spectrum of 7 ^ Vel is composed 
of a constant cooler component (interpreted as arising fur¬ 
ther in the winds), a constant hotter component (which 
would not be sensitive to absorption changes), and a vary¬ 
ing medium-temperature component strongly affected by 
changes in absorbing column. However, a large reduction 
in mass-loss rate was required to fit the observed absorp¬ 
tion value (Schild et ah, 2004). 

The lightcurve of V 444 Cyg also displays a spectacular 
modulation in the soft band (0.4-2. keV), with a peak oc¬ 
curring after the 0 -star conjonction and a decline steeper 
than the rise (Lomax et ah, 2015). This peak can be ex¬ 
plained by three factors: the smaller absorption by the WR 
wind as the line-of-sight avoids the densest parts of the WR 
wind as the shock cone wrapped around the 0 -star turns 
towards the observer, the lower absorption by the more 
tenuous 0 -star wind (filling the shock cone) entering into 
the line-of-sight, and the addition of the 0 -star’s intrinsic 
X-ray emission. Note that a simple model considering only 
absorption by the WR wind is already able to explain the 
presence of a peak in the soft band. The shift of the soft 
emission peak with respect to conjunction time is linked 
to the emitting regions not being exactly centered on the 
binary axis, suggesting aberration of the shock cone by 
Coriolis deflection. The asymmetry of the peak can then 
be reproduced by considering a brighter emission on the 
leading side of the collision zone. Additional evidence for 
the role of Coriolis deflection in this system comes from 
the secondary eclipse of the hard X-ray emission that is 
not centered on the time of optical eclipse(Lomax et ah, 
2015). Finally, the wide shock opening angle inferred from 
the XMM-Newton data indicates that radiative braking 
must be at work in this system (Lomax et ah, 2015), as 
proposed by Cayley et al. (1997). 

In WR 25 (WN 6 ha -H 04f, e ~ 0.5, Porb = 208d), the 
flux at energies below 2 keV rises as the stars move towards 
periastron (separation effect, see below) but this rise is 
abruptly stopped by a sharp decline occurring as the WR 
star passes in front, nearly at the same time as periastron 
(Pandey et ah, 2014). The lightcurve then slowly recovers, 
as the WR wind slowly quits the line-of-sight. 

In WR 140, the situation is similar: a gradual increase 
of the flux as the stars move towards periastron is inter¬ 
rupted by a sharp decrease due to absorption when the 
O star passes behind the WR as seen from Earth (Corco¬ 
ran, 2003; Corcoran et al., 2011). The lightcurves appear 
compatible with that expected from a point source at the 
stagnation point rotating with the binary period (Corco¬ 
ran et ah, 2011 ), but the observed value of absorption re¬ 


mains lower than predicted (Pollock et ah, 2005). Other 
examples of sharp drops in the X-ray flux when the WR 
wind comes into view include WR 21a (WN -|- O, e = 0.7, 
P’orb = 32 d. Cosset and Naze, 2015) and WR 22 (WN7ha 
-I- O, e ~ 0.6, Porb = 80d. Cosset et ah, 2009). In the 
latter case, the observed value of the absorption remains 
much below the predicted value (Cosset et ah, 2009). 

Two examples of absorption effects have also been found 
in the Magellanic Clouds. HD 5980 (LBV -I- WR, e = 0.3, 
Porb = 19.3 d) was the first extragalactic massive binary 
where X-ray emission from colliding winds was securely 
identihed (Naze et ah, 2007). The system displays a clear 
increase in observed flux when the star with the less dense 
wind passes in front (Naze et al., 2007). The triple sys¬ 
tem HD 36 402 (WC4 -|- O -I- O with orbital periods 3d 
and 4.7 yr, though the orbits are not well known) displays 
a hard, bright, and variable X-ray flux compatible with 
the 4.7yr period (Williams et ah, 2013). The system ap¬ 
pears X-ray faint when it is IR bright. This may be linked 
to the periastron passage as it is the time when episodic 
dust formation occurs. The wind interaction region would 
then be deep inside the WR wind and the resulting strong 
absorption would extinguish the X-ray emission. 

Beside WR systems, another example of the effect of 
changing wind attenuation is provided by the non-eclipsing 
binary HDE 228 766 (Porb = 10.7d, e = 0). This sys¬ 
tem hosts an 08 H primary and a secondary star which 
is probably in an intermediate evolutionary stage between 
an Of supergiant and an WN star. Rauw et al. (2014a) 
analysed XMM-Newton spectra at both conjunctions and 
at one quadrature phase to measure the variations of the 
line-of-sight absorption toward the X-ray emission from 
the wind-wind interaction zone and thus infer information 
on the wind densities of both stars. A toy model based 
on the formalism of Canto et al. (1996) was used to set 
constraints on the stellar wind parameters by attempting 
to reproduce the observed variations of the relative fluxes 
and wind optical depths at 1 keV. The variations of the 
X-ray spectra suggest an inclination in the range 54-61° 
and a value of 7 < 0 . 2 . 

Examples of the effect of changing absorption along the 
line-of-sight can also be found among O -I- OB binaries. 
For instance, Arias et al. (2002) showed that the ROSAT- 
PSPC lightcurve of HD 165 052 (06.5 V -k 07.5 V, e = 
0.09, Porb = 2.95 d) displays a double-peaked structure 
with maxima at quadrature phases in agreement with ex¬ 
pectations from theoretical models for binaries with roughly 
equal wind strengths (Pittard and Stevens, 2002). Another 
case is probably the triple system HD 150 136 consisting 
of an 03 V((f*))-3.5 V((f''')) primary and an 05.5-6V((f)) 
secondary on a circular orbit (Porb = 2.67 d) and an 06.5- 
7 V((f)) tertiary on a much wider orbit (Mahy et ah, 2012). 
Skinner et al. (2005) presented a 90 ks Chandra observa¬ 
tion of this system which displays variability on a timescale 
of less than one day. When folded with the ephemerides 
of Mahy et al. (2012), the lightcurve yields a minimum at 
phase 0 . 0 , i.e. at the conjunction phase with the primary 
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component in front. 

However, the detection of such a modulation is not al¬ 
ways a proof for the presence of additional X-ray emis¬ 
sion coming from the interacting wind region. Indeed, in 
the case of the contact binary 29 CMa (08Iaf -I- 09.51, 
e = 0, Porb = 4.39 d) the POS'HT-PSPC lightcurve shows 
a roughly sinusoidal modulation with a maximum when 
the secondary star is in front and a minimum when the 
primary is in front. Berghoefer and Schmitt (1995) sug¬ 
gested to explain this result through the sole effect of the 
wind opacities on the intrinsic emission of each compoent 
without the need to invoke additional X-ray emission from 
a wind-wind interaction. 

There is also at least one counterexample where a mod¬ 
ulation of the X-ray emission due to a changing line-of- 
sight optical depth is expected, but is not observed: In 
WR 79 (WC6-7 -h 05-8, e = 0, Porb = 8.9 d), the X-ray 
emission in an XMM-Newton monitoring was found to be 
constant despite strong expected changes in this asymmet¬ 
ric system (Cosset et ah, 2011). High spatial resolution 
Chandra data revealed the presence of two nearby sources, 
within the XMM-Newton PSF hence affecting the XMM- 
Newton results (Cosset et ah, 2011; Zhekov, 2012). How¬ 
ever, while any variability analysis should be considered as 
preliminary because of that contamination, it would be an 
extraordinary coincidence that a change in WR 79 would 
exactly be compensated by a change in those sources in 
the XMM-Newton monitoring: the observed constancy is 
thus more secure, but is difficult to interpret in the context 
of colliding winds. 

Provided the orbital inclination is right, one could ex¬ 
pect eclipses of the wind interaction zone. Indeed, the 
strongest, head-on wind collision occurs at the so-called 
stagnation point, located on or near the binary axis. This 
central region is regularly occulted by the stellar bodies 
in high-inclination systems. This is the case of V 444 Cyg 
where monitoring with XMM-Newton revealed the pres¬ 
ence of eclipses in the hard (> 2keV) band (Lomax et ah, 
2015, see above). 

In WR 20a (WN6ha -h WN6ha, e = 0, Porb = 3.7d), 
the two stars are identical hence no absorption change is 
expected but the system undergoes photometric eclipses 
because of its high inclination (about 75°, Rauw et ah, 
2005, and references therein) hence one could also expect 
X-ray eclipses. However, a set of Chandra observations 
showed no eclipses in the X-ray domain, on the contrary, 
the system appeared brighter at conjunction phases than 
at phases just before (Naze et ah, 2008). The absence of 
eclipses can only be understood if the X-ray emission zone 
is large while the fainter emission near quadrature could be 
explained by (self-)absorption in the shocked winds (Naze 
et al., 2008). Hydrodynamic models of the system by 
Montes et al. (2013) confirmed these conclusions. Because 
of the small separation between the stars, the winds collide 
while they are still accelerating. Along the binary axis, the 


winds thus collide at low speed and the innermost part of 
the wind interaction region is most probably highly radia¬ 
tive (Montes et al., 2013). The plasma there is thus likely 
too cold to contribute significantly to the X-ray emission. 
The main contribution to the X-ray emission thus arises 
from the outer parts of the wind interaction zone and this 
could explain the absence of eclipses in X-rays as found 
by Naze et al. (2008). The models of Montes et al. (2013) 
also predict sharp decreases in the X-ray fluxes at quadra¬ 
ture phases because of the large column density along the 
contact discontinuity. However, these simulations neglect 
the effect of Coriolis forces which bend the interaction re¬ 
gion hence smearing out the events predicted by Montes 
et al. (2013) and leading to a rather broad maximum of 
the column density (Lemaster et al., 2007). 

In the eclipsing binary CPD-41° 7742 (09 V -h Bl- 
1.5 V, e = 0.02, Porb = 2.44 d) XMM-Newton EPIC data 
reveal a peculiar lightcurve. The X-ray flux is constant at 
a low level over half of the orbit when the rear side of the 
B-star is seen. The X-ray flux then increases except for a 
clear drop that matches the optical secondary eclipse (Sana 
et ah, 2005). This lightcurve is interpreted as the signa¬ 
ture of the primary wind crashing onto the photosphere of 
the secondary or, in case radiative inhibition and radiative 
braking are efficient, a wind interaction zone located very 
close to the surface of the secondary. This situation gen¬ 
erates a hot plasma close to the photosphere of the B-star 
and an effect similar to lunar phases explains the observed 
gradual rise and decline in flux (see Sana et al., 2005). 

3.2.2. Orbital variations due to changing separation 

In addition to the changing line-of-sight column densi¬ 
ties, eccentric binaries can also display changes of the in¬ 
trinsic colliding wind X-ray emission because of the chang¬ 
ing separation d between the stars. For fully adiabatic 
wind interactions, the separation between the stars is usu¬ 
ally high and the winds should thus collide with velocities 
close to Voo at all orbital phases. Therefore, Equation (3) 
predicts variations of the X-ray flux as 1/d. This results 
in a brightening as the stars move towards periastron, and 
a decrease of the X-ray flux as the stars separate again. 
This effect should be most prominent at higher energies, 
as these photons are less affected by absorption. In sys¬ 
tems with purely radiative wind interaction zones, the in¬ 
teraction occurs at lower wind speeds around periastron, 
therefore leading to a softer and also fainter emission. 

Probably one of the best examples of an adiabatic wind- 
wind interaction displaying a l/d variation is provided by 
the non-thermal radio emitter CygOB2 (see Fig. 2). 
Naze et al. (2012b) present XMM-Newton and Swift obser¬ 
vations collected around the 2011 periastron passage. The 
only slight departure from the 1/d relation concerns the 
hard X-ray flux which falls slightly short of the expected 
value at periastron. Naze et al. (2012b) suggest that this 
may be due to the collision becoming slightly radiative as 
a result of the increased wind density at periastron. Naze 
et al. (2012b) also found that the emission was somewhat 
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Figure 2: Variation of the XMM-Newton EPIC-pn spectrum of 
CygOB2 #9 between 0 = 0.190 (black data points) and 0 = 0.997 
(red data points). Note the strong variation in flux level between 
the two observations, the one near periastron being about a factor 
4 higher in good agreement with the 1/d variation between both 
observations (Naze et al., 2012b). 

softer near periastron and suggested that this might hint at 
the effects of radiative inhibition or braking (Stevens and 
Pollock, 1994; Gayley et ah, 1997) becoming more efficient 
as the separation between the stars decreases. This asser¬ 
tion was confirmed by detailed 3-D adaptive mesh refine¬ 
ment hydrodynamic simulations of Cyg OB2 #9 account¬ 
ing for wind acceleration, radiative cooling, and the orbital 
motion of the stars presented by Parkin et al. (2014). In 
these simulations, wind acceleration was found to be in¬ 
hibited at all phases by the radiation field of the opposing 
star. Naze et al. (2012b) showed that the X-ray emission 
of this system is intrinsically quite hard {kT > 1.5 keV) as 
expected for an adiabatic wind interaction zone. This hot 
shocked plasma of the wind interaction zone was also found 
to contribute significantly to the free-free radio emission of 
the system (Blomme et ah, 2013). 

The non-thermal radio emitter 9 Sgr is another exam¬ 
ple of a highly adiabatic wind interaction zone. However, 
Rauw et al. (2015) found that the hard (2.0 - 10.0 keV) 
X-ray flux of 9 Sgr at periastron is lower than expected 
from a simple 1/d scaling relation. Since the stars have 
a weaker radiation field and are much farther apart from 
each other than in CygOB2 #9, Rauw et al. (2015) argue 
that radiative inhibition or braking are unlikely to play a 
role. They suggest instead that shock modification by rel¬ 
ativistic ions (Pittard and Dougherty, 2006) might be at 
work in this system. 

Other systems that show hints for a 1/d variation of 
their X-ray flux, at least at high energies, are HD 93403 
(05.51 -t- 07V, e = 0.23, Porb = 15.1 d, Rauw et ah, 
2002), HD 93 205 (03V -b 08V, e = 0.37, Porb = 6.1 d, 
Antokhin et ah, 2003) and WR 25 (Pandey et ah, 2014). In 
HD 93 403 and WR 25, the variations in the softer energy 
bands are contaminated by absorption effects. Moreover, 
at least in the case of HD 93 403, the offset due to the in¬ 
trinsic emission of the stars yields changes of the X-ray flux 


that are lower than expected from a simple 1/d scaling. 

However, the predicted 1/d variations are not always 
observed. In WR 22 and 7 ^ Vel, the hard part of the 
spectrum remains constant despite the large variation in 
separation due to a high eccentricity (Gosset et al., 2009; 
Rauw et ah, 2000; Schild et ah, 2004). Parkin and Gosset 
(2011) presented tailored hydrodynamic models of WR 22. 
Both instantaneously accelerated and radiatively driven 
stellar winds were considered and compared against the 
observations of Gosset et al. (2009). Neither of the models 
actually reproduced the observed lightcurve well. 

In the same context, repeated monitoring campaigns of 
WR 140 prooved the existence of variations of its intrinsic 
flux. However, these variations by about a factor three 
remain well below those expected (a factor 19) from the 
1/d relation for such a large eccentricity (Pollock et ah, 
1995; Zhekov and Skinner, 2000; Gorcoran et ah, 2011; 
Sugawara et ah, 2011). As pointed out above, Pittard and 
Dougherty (2006) suggested that, in the case of WR 140, 
this situation results from the impact of particle accelera¬ 
tion on the properties of the shocks. 

For shorter period eccentric systems, there are strong 
deviations from the 1 /d relation and the wind interaction 
zone can switch from mostly adiabatic near apastron to 
highly radiative near periastron. As a result, theoretical 
models of Pittard and Parkin (2010) predict that the X-ray 
flux at a given orbital phase depends on the history of the 
plasma heated at earlier orbital phases. There are now 
a number of systems where these predictions have been 
verified by observations. To identify such a behaviour re¬ 
quires a rather good phase coverage of the orbital cycle. 
For instance, Gazorla et al. (2014) analysed a sample of 
XMM-Newton and Swift observations of the core of the 
Gygnus OB association, significantly improving the phase 
coverage of the 21.9 d eccentric 06If -b 05.5HI(f) binary 
GygOB2 # 8 a (e = 0.21). The X-ray emission clearly dis¬ 
plays phase-locked variations that are anticorrelated with 
the changes of the non-thermal radio emission of this sys¬ 
tem. In a system such as CygOB2 # 8 a, where the wind 
interaction zone is at least partially radiative, the collision 
occurs at higher speeds at apastron, when the winds have 
enough room to accelerate. This leads to a stronger shock 
at apastron, and thus a harder X-ray emission. However, 
the X-ray flux does not peak at apastron, but at an inter¬ 
mediate phase of ^ ^ 0.8. A plot of the X-ray flux ver¬ 
sus the orbital separation (Fig. 3 of Gazorla et ah, 2014) 
revealed striking similarities with Fig. 19 of Pittard and 
Parkin (2010). The predicted hysteresis behavior is clearly 
seen with the emission being harder at phases when the 
stars get closer than when they separate again. 

Similar effects are present in HD 152 248 (07.5(f) HI 
-b 07(f) HI, e = 0.13, Porb = 5.816 d) and HD 152 218 
(09IV -b 09.7 V, e = 0.26, Porb = 5.60 d) as illustrated in 
Fig. 3. Sana et al. (2004) found the X-ray flux of the highly 
radiative system HD 152 248 to display a clear, asymmetric 
modulation with orbital phases, reaching a maximum at 
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(j) = 0.66. In the case of HD 152 218, Sana et al. (2008) 
also found that the X-ray flux and its hardness increase 
near apastron. 

Yet another example of an hysteresis is found in the X- 
ray lightcurve of WR 21a (WN -|- O, e = 0.7, Porb = 32 d). 
In this system an increase of the flux is detected towards 
periastron, but the flux variations are clearly dependent on 
both separation and phase (Gosset and Naze, 2015). Fur¬ 
thermore, just before periastron, the X-ray flux suddenly 
decreases. This change in behaviour cannot be due to the 
sole effect of an increasing absorption as it is also seen at 
hard energies which are far less affected by photoelectric 
absorption. This situation instead suggests a crash of the 
wind interaction zone onto the surface of the 0-star (Gos¬ 
set and Naze, 2015). 

There are a number of systems where X-ray observa¬ 
tions suggested the presence of wind interactions, but ac¬ 
curate orbital elements are so far lacking. This is the case 
of the triple or quadruple system CygOB2 (Cazorla 
et al., 2014). The variations of the X-ray flux indicate 
that they do not arise from the short-period ( 6.6 days) 
06.5-71 -I- OB-Ofpe/WN9 inner binary, but rather from 
the collision between the winds of the inner binary and 
that of the tertiary component orbiting it with a 6.7yr 
period as derived from the modulation of the non-thermal 
radio emission (Kennedy et al., 2010). Yet, there are cur¬ 
rently huge uncertainties on the orbital elements of the ter¬ 
tiary which prevent a more thorough interpretation (Ca¬ 
zorla et al., 2014). The same is true for WR 48a (Zhekov 
et al., 2014a,b) and WR 65 (Oskinova and Hamann, 2008): 
variations of the X-ray emission have been observed, but 
are difficult to interpret without a detailed knowledge of 
the orbit. 

In other systems, the scarce number of data prevents a 
detailed investigation. For instance, XMM-Newton data of 
WR 137 suggest some variations compared to the Einstein 
or ROSAT observations though all observations were taken 
at similar separations (Zhekov, 2015). 

3.3. Observational insight into shock physics 

As pointed out in Sect. 2 . 2 , in wide systems, the shock- 
heated plasma could be out of ionization equilibrium and 
the electron and ion temperatures are likely non-equal. 
For WR 140, Zhekov and Skinner (2000) found that mod¬ 
els with non-equal electron and ion temperatures were fa¬ 
vored. This finding was confirmed by Pollock et al. (2005) 
on the basis of theoretical considerations (at apastron, the 
Coulomb interaction length is similar to the orbital sepa¬ 
ration) , but also from the fact that line width seems to in¬ 
crease with ionization potential. This latter result suggests 
that Ne and Mg lines are produced more centrally (closer 
to the binary axis where the spread in velocities is expected 
to be smaller) than the Si or S lines. In ionization equilib¬ 
rium, one would rather expect the hottest plasma (hence 
emission lines associated with ions of higher ionization po¬ 
tential) to arise closer to the binary axis. A similar prob¬ 


lem is present when comparing the width of He-like triplets 
to lines of H-like ions, again suggesting that they do not 
arise from the same region. Pollock et al. (2005) thus con¬ 
cluded that ionization balance is not reached. These au¬ 
thors also reported the simultaneous presence of lines from 
many different ionization stages of iron (Fexvii - Fexxv) 
in the ChandraAAEAG spectra between 10.6 and 15.0 A. 
Pollock et al. (2005) suggested that there is no efhcient 
electron heating of the pre-shock gas because of the mag¬ 
netic field at the shock (WR 140 is a non-thermal radio 
emitter) but no detailed calculations were presented. 

In a similar manner, Zhekov (2007) suggested that, in 
view of timescales, ionization equillibrium should not be 
reached in the very wide system WR 147 (WN 8 -I- OB). He 
noted that non-equilibrium ionization, collisionless shock 
models give a better fit to the data. 

For 7 ^ Vel, the scaling relations indicate that the shocked 
plasma should be in ionization equilibrium. However, the 
analysis of Chandra data by Henley et al. (2005) indicated 
that Mg lines are formed in hotter plasma located closer to 
the 0-star than Si lines. This is a typical signature of non¬ 
equilibrium ionization. This result was however not recov¬ 
ered in the analysis of Skinner et al. (2001) of the same 
data, as the authors used different calibration files and 
background subtraction: the presence of non-equilibrium 
ionization in 7 ^ Vel remains thus currently an open issue. 

A very peculiar situation concerns the wind interac¬ 
tion in Plaskett’s Star (HD 47129, 061 -|- 07.51, e = 0, 
T^orb = 14.4 d). Combining XMM-Newton observations 
with ROSAT data, Linder et al. (2006) found a modulation 
of the X-ray emission with orbital phase with a minimum 
when the primary star is in front, although there was quite 
substantial scatter in the lightcurve. This suggests quite 
a strong occultation effect by the body of the primary 
star in this non-eclipsing binary and thus implies that the 
wind interaction zone is located closer to the surface of the 
primary than one would expect from the wind momentum 
balance. Optical spectra of HD 47129 presented by Linder 
et al. (2008) revealed an annulus of circumstellar material, 
likely due to a confinement of the wind of the secondary. 
The primary wind would then collide with this confined 
wind, leading to strong X-ray emission at the contact be¬ 
tween the confined secondary wind and the incoming pri¬ 
mary wind and shifting the interaction zone closer towards 
the primary than one would expect without confinement of 
the secondary wind. Whilst Linder et al. (2008) attributed 
this confinement to the fast rotation of the secondary, sub¬ 
sequent spectropolarimetric observations by Grunhut et al. 
(2013) led to the discovery of a magnetic field of dipole po¬ 
lar strength 2.85 ± 0.5 kG, suggesting that this magnetic 
field is responsible for the confinement of the secondary’s 
wind. Therefore, Plaskett’s Star likely features a complex 
combination of X-ray emission from the shocks within the 
secondary’s wind as it is confined by the magnetic field as 
well as from the collision between the primary wind and 
the magnetosphere of the secondary. To date, this is the 
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Figure 3: Left: total X-ray flux of HD 152 248 observed with XMM-Newton as a function of orbital separation (Sana et al., 2004). Right: 
same for the medium-band flux of HD 152 218 (Sana et al., 2008). The hysteresis-like behaviour is clearly seen in both cases. The orbital 
phases of the observations are indicated by the labels. 


only system known to exhibit such a configuration. 

3.4- High resolution X-ray spectroscopy 

Chandra and XMM-Newton have initiated the era of 
high-resolution X-ray spectroscopy for moderately bright 
X-ray sources. For massive stars, HETG and RGS spectra 
unveiled for the first time some details of the morphology 
of the spectral lines (see e.g. Giidel and Naze, 2010, for a 
review). 

Only few interacting wind systems have been observed 
at high spectral resolution in X-rays. This is mostly due 
to the limited sensitivity of current facilities, but clearly 
this is a promising route for future observations. Indeed, 
important information can be derived from high-resolution 
spectra, including the chemical composition of the emit¬ 
ting plasma, the position of the emission region (through 
the analysis of He-like triplets), and the velocity of the 
emitting plasma (through line width and shift). For ex¬ 
ample, in single massive stars, the forbidden lines of He- 
like triplets are suppressed through UV pumping of the 
upper level of the transition as the X-ray emission occurs 
in a region close to the star, hence flooded by intense UV 
radiation. Wind-wind collisions occur farther away from 
the photospheres and the forbidden lines therefore remain 
strong (e.g. WR 147, Zhekov and Park 2010a; WR 48a, 
Zhekov et al. 2014a; WR140, Pollock et al. 2005; 6 Muscae, 
Sugawara et al. 2008; y^Vel, Skinner et al. 2001; Schild 
et al. 2004) , which is a strong argument in favour of X-ray 
emission arising inside the colliding wind region. 

Furthermore, while detailed line profile analysis is not 
yet possible, the line shift and width have been recorded in 


a few cases. Widths are generally in the 1000-2000 km s“^ 
range (Pollock et ah, 2005; Sugawara et al., 2008; Zhekov 
and Park, 2010a; Zhekov et ah, 2014a), i.e. of the same or¬ 
der of magnitude as wind speeds, as expected. Line shifts 
are directly related to the orientation of the colliding wind 
zone. For example, in the triple system 0 Mus (WC5-6 
-I- 06-7V -I- OBI, Porb = 19d), lines appear redshifted 
(Sugawara et al., 2008): if the orbital solution of the close 
binary is correct, then these redshifted lines cannot be 
associated with a shock cone within the binary, folded 
around its O component because its orientation at that 
time would rather yield blueshifted lines. Instead, these 
redshifts would arise in a collision between the combined 
wind of the binary and the wind of its optical compan¬ 
ion, an 0-type supergiant (Sugawara et ah, 2008). Still, 
the ephemeris must be checked and the line shifts must 
be monitored closely before definitive conclusions can be 
drawn in this case. 

A very special situation is provided by the very wide 
system WR 147. To date, this is the only example where 
the spatial extension and location of the X-ray source 
were measured using high-spatial resolution Chandra data. 
First evidence for spatially extended emission associated 
to a wind interaction zone was reported by Pittard et al. 
(2002). Subsequent deeper observations by Zhekov and 
Park (2010b) revealed further details: a northern source 
which is extended and coincident with the bowshock ob¬ 
served at radio wavelengths, and a southern one that co¬ 
incides with the position of the WN8 star. In view of its 
location and spectrum, the former component was there¬ 
fore considered as the signature of the wind-wind colli¬ 
sion between the WR and the OB companion. However, 
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the latter component displays a bright, hard, and vari¬ 
able emission untypical of late WN stars: Zhekov and 
Park (2010a) therefore proposed that the WR actually is 
a close binary where a second wind-wind collision takes 
place. In WR 147, the lines from the northern compo¬ 
nent appear blueshifted, indicating a colliding wind zone 
located between the observer and the WR (Zhekov and 
Park, 2010a). However, owing to its very long period, the 
orbit of WR 147 is essentially unknown. 

The orbit of WR 140 is far better determined and the 
two high-resolution spectra presented by Pollock et al. 
(2005) reveal line variations clearly associated to the ge¬ 
ometry of the shock zone. Before periastron, the shock 
cone is opened towards us and the lines naturally appear 
blueshifted; after periastron, the cone has turned away 
and the lines are thus slightly redshifted (and somewhat 
broader). 

In the poorly-known episodic dust maker WR 48a, spec¬ 
tral lines seen in a high-resolution Chandra spectrum are 
also blueshifted but this result is more difficult to inter¬ 
pret because of the absence of accurate orbital elements 
(Zhekov et al., 2014a). In fact, the Chandra observation 
indicates a larger absorption than what is measured on an 
XMM-Newton observation of WR 48a (Zhekov et al., 2011, 
2014a), suggesting the WR-star was in front of its compan¬ 
ion at the time of the Chandra observation. The observed 
blueshift suggests that the cone is wrapped around the WR 
star (or the star in the system with the wind producing 
the largest absorption). This implies that the companion 
should have a more energetic wind, therefore suggesting 
that WR 48a is a WR -I- WR system, instead of a classical 
WR -I- O binary (Zhekov et al., 2014a). This assumption 
was confirmed by subsequent optical and near-IR observa¬ 
tions that revealed a WC 8 -I- WN 8 h system (Zhekov et al., 
2014b). 

In 7 ^ Vel, no line shift was detected on a Chandra grat¬ 
ing observation taken near X-ray maximum (Skinner et al., 
2001; Henley et al., 2005), which is at odds with expec¬ 
tations. At that phase, one expects the shock cone to 
be opened towards Earth, which is demonstrated by the 
higher apparent flux linked to the lower absorption of the 
0-star wind (see above). To get a negligible line shift un¬ 
der these conditions, the only possibility is to have a wide 
opening of the shock cone (half opening angle > 80° Hen¬ 
ley et al., 2005). This could be due to radiative braking 
by the 0 -star’s radiation, but a wide cone is incompati¬ 
ble with the X-ray lightcurve: the short duration of the 
low absorption episode, occurring when the line-of-sight 
crosses the wind of the 0 -star clearly requires a narrow 
opening (Willis et al., 1995; Henley et al., 2005, half open¬ 
ing angle ^25°). A possible solution to the problem could 
be contamination of the spectral lines by the intrinsic X- 
ray emission of the 0 -star which should be essentially un¬ 
shifted. 

The X-ray spectra of several carbon-rich Wolf-Rayet 
binaries exhibit radiative recombination continua (RRCs) 


mainly associated with carbon. This is the case for 7 ^ Vel 
(Schild et al., 2004) and 9 Mus (Sugawara et al., 2008)®. 
For the latter system, this feature appears redshifted. The 
presence of RRCs requires the presence of low-temperature 
gas that is ionized by the X-ray emission (probably from 
the colliding wind interaction). This cool gas could be 
associated with material downstream of the shock cone. 

4. Conclusions and future perspectives 

Over the past fifteen years, there has been tremendous 
progress in our understanding of wind-wind interactions in 
early-type binaries. Theoretical models with an increasing 
degree of sophistication have been developed and observa¬ 
tions with the current generation of X-ray observatories 
allowed to study the phenomenon with unprecedented de¬ 
tails over a wide range of the parameter space. In par¬ 
ticular monitoring of interacting wind systems with low- 
resolution CCD spectroscopy has revealed the complexity 
of the X-ray lightcurves as well as their energy dependence. 
Some of the observational results, such as the 1/d flux 
modulation in some eccentric adiabatic systems and the 
hysteresis in shorter-period eccentric systems, provide nice 
confirmation of the theoretical predictions. Other findings, 
such as the lack of 1/d variations in some wide eccentric 
systems or the lower than expected level of the colliding 
wind X-ray emission, still challenge the theory. An impor¬ 
tant question to be addressed in future studies concerns 
a better understanding of the ionization conditions (equi¬ 
librium versus non-equilibrium) of the post-shock plasma. 
Another aspect to be investigated is why some massive bi¬ 
nary systems are X-ray dim whereas others are prominent 
X-ray emitters. To address the latter question one needs 
to consider binary systems over as wide a range as possible 
of the relevant parameters (orbital separations, eccentric¬ 
ities, mass-loss rates, wind velocities,...). In addition to 
X-ray monitoring campaigns, such studies also require a 
major observational effort in the optical domain to estab¬ 
lish accurate orbital parameters. 

The coming generations of X-ray observatories, espe¬ 
cially the JAXA-led mission Astro-H (Takahashi et al., 
2014) and the ESA-led observatory Athena (Nandra et al., 
2013), will enable even finer observational studies. In fact, 
these observatories will carry microcalorimeter spectro¬ 
graphs, called the Soft X-ray Spectrometer (SXS, Mitsuda 
et al., 2014) for Astro-Hand the X-ray Integral Field Unit 
(X-IFU, Barret et al., 2013; Raveraet al., 2014) for Athena. 
These instruments will have unprecedented spectral reso¬ 
lution around the 6.7keV Fe xxv K complex, allowing for 
the first time to resolve the line profiles and follow their 
changes. For early-type stars, the Fe K line is only de¬ 
tected in the spectra of colliding wind binaries (e.g. Fig. 2) 


®In addition, hints of an RRC in a lower-resolution Suzaku CCD 
spectrum of WR 140 were reported by Sugawara et al. (2011), al¬ 
though such detections are less secure than those with high-resolution 
spectrographs. 
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or of stars with magnetically confined stellar winds. It is 
thus an efficient tracer of the presence of very hot plasma, 
which - in the case of wind-wind collisions - arises in the 
innermost part of the interaction zone. Phase-resolved ob¬ 
servations of the Fe K lines for a sample of massive bina¬ 
ries (with circular or eccentric orbits) will hence provide 
new insight into the physics of the phenomenon (Sciortino 
et ah, 2013; Rauw et ah, 2016). AthenawiW also allow us to 
enlarge the sample of extragalactic colliding wind binaries 
that can be studied in X-rays. This is especially interesting 
in terms of the theoretically predicted dependence of the 
mass-loss rates on metallicity. Such a dependence could 
affect both the intrinsic level of the X-ray emission of the 
wind-wind interaction zone and the level of photoelectric 
absorption by the unshocked winds. 
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